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Abstract 

We have investigated the impact of long-range forces induced by unparticle operators of scalar, 
vector and tensor nature coupled to fermions in the interpretation of solar neutrinos and KamLAND 
data. If the unparticle couplings to the neutrinos are mildly non-universal, such long-range forces 
will not factorize out in the neutrino flavour evolution. As a consequence large deviations from the 
observed standard matter-induced oscillation pattern for solar neutrinos would be generated. In 
this case, severe limits can be set on the infrared fix point scale, A u , and the new physics scale, M, 
as a function of the ultraviolet (dijv) an d anomalous (d) dimension of the unparticle operator. For a 
scalar unparticle, for instance, assuming the non-universality of the lepton couplings to unparticles 
to be of the order of a few per mil we find that, for d\jy = 3 and d = 1.1, M is constrained to 
be M > O(10 9 ) TeV (M > O(10 10 ) TeV) if A u = 1 TeV (10 TeV). For given values of A u and d, 
the corresponding bounds on M for vector (tensor) unparticles are ~ 100 (~ 3/y A„/TeV) times 
those for the scalar case. Conversely, these results can be translated into severe constraints on 
universality violation of the fermion couplings to unparticle operators with scales which can be 
accessible at future colliders. 
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I. INTRODUCTION 



We are accustomed to describe fundamental interactions in nature in terms of quantum 
fields which manifest themselves as particles. Our most successful description of all the 
available experimental data today is given by the Standard Model (SM) which is based on 
such quantum fields. Extensions of the Standard Model, such as supersymmetric, grand- 
unified, or extradimensional models, also rely on these ideas and give rise to a plethora of 
new particles that may be discovered soon at the LHC. However, other types of quantum 
fields may be present in Nature. 

Recently, Georgi [l|, [2| proposed a novel scenario where no new particle states emerge 
but where scale invariant fields with non-trivial infrared (IR) fixed point can be part of 
physics above the TeV scale. If these fields couple to SM operators they can give rise to 
effective interactions at lower energies that may have impact on particle phenomenology. He 
denominated this scale invariant matter coupled to SM fields at lower energies unparticle, 
as it is not constrained by a dispersion relation like normal particles. 

Since then there has been a number of papers in the literature devoted to map the 
unusual observable consequences of unparticles to experiments. QED bounds to unparticle 
interactions have been obtained in [3|. Signals of unparticle have been studied in past and 
future collider experiments (3, 0]. Possible unparticle contributions to CP violation deep 
inelastic scattering 0], lepton flavour violating processes j8| as well as to hadron mixing 



and decay [9|, neutrino interactions |10l. Illl] and nucleon decay [12| have been investigated. 



Several cosmological and astrophysical bounds have also been devised [13|, [14|, [15[ . 

Another striking spin-off of the unparticle idea is the advent of long-range forces operating 
at macroscopic distances and generally governed by a non-integral power law. This unique 
behavior was first pointed out in Ref. [16|, in connection to microscopic spin-dependent 
interactions between electrons and has also been investigated in connection to deviations 
of the Newtonian gravitational inverse squared law mediated by a tensor-like ungravity 



interaction [17j. Limits on vector-like unp article interactions coupled to baryon and lepton 
flavour numbers were established in 18| by comparing with results from torsion-balance 
experiments. 

Neutrino oscillations have been shown to be a very powerful tool in constraining new 
dynamics. Being an interference phenomenon it can be sensitive to very feeble interactions 
as their effect enters linearly in the neutrino propagation. Furthermore the effect of new 
fermionic interactions can be enhanced when neutrino travels regions densely populated by 
those fermions. In this way it has been shown that neutrino oscillation data can place strin- 
gent bounds on long-range forces coupled to lepton flavour numbers j26|, 27, 28[. Therefore 
it is expected that they can provide very strong limits also on the couplings of unparticles 
to fermionic currents. 

Assuming that scale invariance is not broken in the IR we derive in this paper, from 
solar and reactor neutrino experimental data, limits on scalar, vector and tensor unparticle 
operators which mediate long-range interactions through their couplings to fermionic flavour- 
conserving currents. In Sec. [Ill we briefly present the formalism of unparticle dynamics 
relevant to the present study. We present the effects of the induced long-range forces in 
solar neutrino oscillations in Sec. IHIl We derive the bounds from the no observation of such 
long-range effects from the global analysis of solar and KamLAND data in Sec. [IV] and in 
the final Section we discuss our conclusions. 
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II. UNPARTICLE FORMALISM IN BRIEF 



We will consider hypothetical fields of a hidden sector of a theory with a non-trivial IR 
fixed point which can interact with SM fields at very high energies by the exchange of a 
particle with a large mass M as proposed in Ref. So the ultraviolet (UV) couplings of 
this hidden sector operators C?uv of dimension c?uv with SM operators Osm of dimension 
<1sm take the form 

(1) 



][{dij V +d SM - 

Scale invariance emerges at a lower energy scale A u . As the hidden sector runs towards 
the IR fixed point, these couplings are so to speak sequestered, their scaling dimension is 
increased by the anomalous dimension d from the hidden sector dynamics. From this point 
on the hidden sector operator matches onto an unparticle operator 0\j whose interactions 
with the SM operators are described by the non-renormalizable effective lagrangian 

^d uv -d 

Cu M d v U v+ d SM -4 °SM Ojj, (2) 

where C u is a dimensionless constant of order 1. 

Our starting point will be the lagrangian which contains the effective flavour diagonal 
interactions which satisfy the SM gauge symmetry, between scalar, vector and tensor un- 
particles and the SM fermions / 

where the unparticle operators are taken to be hermitian. Additional four fermion contact 
interactions could also arise as discussed in Ref. [l9j]. I n the absence of cancellations they 
will add up to the effects discussed here making the unparticle effects larger. In what follows 
we will conservatively neglect the contributions from contact interactions. 
The natural size of the A{ 17 t coefficients is 



duv— 1 

M 



^S,V,T - ^S,T,V ( T7 ) ' ( 4 ) 



where Cg TV are again of order 1. 

Scale invariance determines the propagation properties of these operators up to a normal- 
ization factor. Using the same convention as in Ref. [I, H[ , including the corrections given in 
Tif. we obtain that the propagators corresponding to scalar, vector and tensor unparticle 



operators are correspondingly 



A(P 2 )= A f 7 A -P 2 ) d ~\ (5) 

v ; 2sin(d7r) v ; v ; 

[A(P 2 )U = A(P 2 )^, (6) 

and [A(P 2 )] W = A(P 2 )T^ pa . (7) 

(8) 

P P 

I ^ u lc\\ 

?V = -fiW + a —j5T ' ( 9 ) 
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where a = 1 for transverse Oy and a = 2 ^_ 1 2 - > in conformal field theories 19| 
For transverse and traceless operators 



T, 



/iv, per 



"p,a"up 



while in conformal field theories 



T, 



[ii/, pa 



+4 



I'U 

(d-2) 
d 

(d-2) 



k,,k„ 



d(d 

In the equations above 



/'/-' 



9V 



A- 2 

kpk a 



[4-d(d+l)] 
2d(d- 1) 



A- 2 



2 

T^Tf fiv^ pa 



9pv9 pa 



V 



(10) 



(11) 



(d 2)(d 3) kpk v kpk a 



d(d 



(k 



2\2 



16tT 



T(d- 



(12) 



(13) 



(27r) 2d r(d- l)T(2d) ' 

These massless unparticle operators induce long-range forces between fermions. In par- 
ticular, the effective lagrangian describing the long-range force generated by a static distri- 
bution of fermions / with number density n./(r), assumed to be spherically symmetric and 
of radius R, on a relativistic fermion /' is given by 



C 



C S s C f s CyCy j, Q , + C S T C^ rrifE' g ji^q j, 



in 



Air 



4tt A 2 



W f (r) 



(14) 



where 5 = 2/3 for a transverse and traceless unparticle tensor operator, while in conformal 
field theories B = 1 — ^j^TTp- The potential function, irrespectively of the Lorentz struc- 
ture of the unparticle operator, can be shown to have a universal form as a function of the 
distance r from the center of the matter distribution 



WAr 



2(d L 



-1) 



2 T[d + ^}T[d 



2 3 
x - 
r i r i 

2 



tt 2 ^- 1 ) T[2d] 



1 



(RA U 



2{d-l) 



\x-y 



3 ~ 2d ] dy if d ^ | 



(^+y) 2 

(x-y) 2 



(15) 



if d 



For d — > 1 one recovers the Coulomb-type potential characteristic of a massless particle 
exchange. This potential function is always finite for d < 2 whereas for d > 2 there are 
infinite contributions that appear when x = y. In order to render a finite answer for d > 2 
this infinite contribution has to be canceled out by the introduction of some counterterm. 
In principle, this could be done by means of the effective four fermion contact interactions 
whose couplings gets renormalized. Technically this is equivalent to remove the point x = y 
from the integration in Eq. (1151) . In principle a finite contact contribution to the potential 
can remain. As mentioned above in what follows we will neglect this contribution. In 
practice we will only be quantitatively discussing the results for d < 2, since, as we will see 
below, the effects become very feeble for larger dimensions. 
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If the theory is not only scale invariant but also conformal invariant the dimension d of 



the unparticle operator is bounded from unitarity constraints For scalar, vector and 
symmetric tensor operators d > 1, 3 and 4, respectively. We will relax this constraint in our 
phenomenological analysis. 

One must also notice that the Higgs coupling to scalar unparticles generally breaks scale 



invariance at a scale close to the electroweak scale [22], |23|, |24j, |25[ . Therefore for the scalar 
unparticle effects to be relevant at solar and reactor energies (~ MeV) discussed here, one 
must assume the Higgs-scalar unparticle couplings to be suppressed. For vector and tensor 
unparticles, the Higgs only couples to higher dimension operators which relaxes the bound 
on the scale of scale-invariance breaking. 



III. UNPARTICLE EFFECTS IN NEUTRINO OSCILLATIONS IN THE SUN 

Under the assumptions discussed in the previous section, unparticles mediate a long-range 
force coupled to fermions in the way described by Eq. (JHj). Such force will affect neutrino 
oscillations in matter in a similar fashion to the leptonic forces discussed in Ref. 



26 



In our description of the effect of unparticles in the evolution of neutrinos in matter we 
have to take into account that conventional mass-induced neutrino oscillations describe very 

~M, M, M, M, and reactoril, M 



well atmospheric [3j, K2K [38(, MINOS |39|, solar 
neutrino data as long as all three neutrino flavours must participate of oscillations [41[ . The 
mixing among neutrino mass eigenstates and flavour eigenstates is encoded in the 3x3 
lepton mixing matrix that in the standard parametrization takes the form 



U 




c 2 i S12 
-S12 C12 
1 



(16) 



where c,- 



cos 6ij and Sjj 



sin 6ij . 



The neutrino oscillation data is consistent with the following hierarchy of the neutrino 
mass squared differences 



Am 







Am 21 <C I Am 



311 



\Am 2 32 \ 



Am 



Recent KamLAND |35l and MINOS |3J 



results lead to Am\ xl 



atm ' 



\Am\ x \ 



(17) 

0.03 HH, so for 



KamLAND and solar neutrinos, the atmospheric neutrino oscillation scale contributions are 
completely averaged and the interpretation of these data in the three neutrino oscillation 
framework depends basically on Am^, #12 and #13. In contrast, atmospheric and long 
baseline neutrino oscillation experiments are driven mostly by Am^, #23 an d #13- Moreover, 
the non observation of oscillations in the CHOOZ reactor experiment [40j imply that #13, the 
mixing angle connecting solar and atmospheric neutrino oscillation channels must be very 
small, i.e., sin 2 # 13 < 0.041 at 3a [ill ]. These considerations make that the 3-v oscillations 
effectively factorize into 2-v oscillations of the two different subsystems: solar plus reactor, 
and atmospheric plus long baseline. 

Because the new couplings are flavour diagonal, their effects in the evolution of atmo- 



spheric neutrinos do not change the hierarchy (j!7p [281 ] . Thus the 1v oscillation factorization 
still holds and as long as we neglect the small #13 we can perform the analysis of solar and 
reactor neutrinos in an effective two-neutrino mixing scenario. In this framework, one can 



5 



easily read from Eq. (fT4j) the effect of the long-range unparticle force in the evolution of 
solar neutrinos as: 



dr 



2E„ 



mi 
m 2 



U 



M UNP (r) 




+ 



(V cc (r 



-VunpW 




(18) 



where E v is the neutrino energy, m 12 are the neutrino masses in vacuum, and U# 12 the mixing 
matrix between neutrino flavour and the vacuum mass eigenstates and v a = 023^ + s 23 z/ T . 
Vcc( r ) = V / 2G j pn e (r) is the Mikheyev-Smirnov-Wolfenstein (MSW) matter potential. 

In Eq. (|T8|) Munp and Vunp(^) depend on the Lorentz structure of the unparticle operator. 
The function W will enter either as an extra mass term (scalar unparticle) or as an addition 
to the MSW potential (vector or tensor unparticle). For scalar unparticles, 



M UNP = M s (r) = 
while for vector unparticles, 

Vu NP (r) = Vy(r) 
and for tensor unparticle operators, 



4ir 



W(r) ~ a s W(r) 



VuNp(r) = 



cuc^ - C 



V • 



4tt 



W(r) ~ a v W(r) 



M 



UNP 



(19) 



(20) 



UNP 



V) = — Vr(r) 



An 



Cp)m e E„ Bw ^ 



m P E v 



A?, 



-olt- 



A 2 , 



W{r), M UNP =((2;) 



where C u s a VT 



23 



^sVt + s \?P V svt- I n a ^ cases, W(r) is given in Eq. (JT51) with 



rif(r) = n e (r) 



C s,v,t „ Cg.y.r ^n(r) 



(22) 



with i? P)n = 1 for scalar and vector unparticles and R P}U = m p ^ n /m e for tensor unparticles. 
Yf(r) is the relative number density of fermion /, it is different for neutrons than for protons 
and electrons because of the change in composition along the Sun radius. Because of this 
composition change, not only the normalization but also the W (r) profile (in r) has some 
dependence on the unparticle couplings. We have neglected this small unparticle effect in 
the profile of the potential and absorbed it in the definition of effective r average couplings: 



oisy,T 



c 'S,V,T\^ 'S,V,T 



r<Va \ 

^S,V,TJ 



-in 



C s,v,t „ Cs,v,t , Y n (r) 

~ K P + ~ e \~U~7Z\' n 



c 



S,V,T 



(23) 



where (yrpj) stands for the average of the relative number densities along the neutrino 
trajectory. 

From Eq. (|2"31 we see that, as expected, the unparticle effect on the solar neutrino evolu- 
tion does not factorize out as long as the unparticle coupling to neutrinos are non-universal. 
As a matter of fact, under the same assumption, solar neutrinos may also invisibly decay 
in their way from the Sun to the Earth as described in Ref. ll|. In what follows we will 
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FIG. 1: Potential function W(r) due to the density of protons (or electrons) in the Sun as a 
function of the distance from the solar center in units of the solar radius, Rq, for various values 
of the dimension d and the mass scale M. We take A u = 1 TeV and d\jy = 3. For comparison, 
we show the value of the W(r) function for a Coulomb-type (d = 1) long-range force with range 
A = 0.1 Rq and with the same strength at the center of the Sun. 



neglect the decay and concentrate on the effects due to the long-range forces, which as we 
will see, turn out to be more constraining. 

In Fig. [1] we show the function W(r) in the Sun as function of the distance from the 
center in units of Rq. As seen in the figure, the dimension d acts as a kind of range of the 
potential. Inside the Sun and up to 0.1 R & , W is constant, independently of d. We notice 
that for any value of d < 1.7, there is always a value of M for which, at any r, the W(r) 
function is always larger than the corresponding one for a d — 1 long-range force of range 
A = 0.1. 

We show in Fig. [2] the survival probability of solar v e at the sunny face of the Earth as 
a function of the neutrino energy E v for some values of the parameters. The probability is 
obtained by numerically solving Eq. (I18p along the neutrino trajectory from its production 
point in the Sun to its detection point on the Earth. We have verified that for the range 
of parameters of interest the evolution in the Sun and from the Sun to the Earth is always 
adiabatic. From the figure we can infer the expected order of magnitude of the bounds that 
can be derived from the solar and KamLAND analysis which we describe next. For the 



scalar case we have chosen mi = so ni2 = \/ Am^. As shown in Ref. 26], the effect of 
scalar interactions becomes stronger as mi grows, so we chose m-i = which gives the most 
conservative constraints. 
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FIG. 2: Survival probability of v e in the Sun as a function of the neutrino energy E v for a scalar 
(upper panel), vector (middle panel) and tensor (lower panel) unparticle force, for various values 



of the parameters. For all curves we have used tan 2 #12 



0.44 and Am^ = 7.9 x 10" 5 eV^. 



IV. CONSTRAINTS FROM SOLAR AND REACTOR NEUTRINO DATA 

We present in this section the results of the global analysis of solar and KamLAND data 
for the unparticle long-range forces discussed in the previous section. Details of our solar 



neutrino analysis have been described in previous papers |41l . |43| . We use the solar fluxes 



from Bahcall, Serenelli and Basu 2005 42J1. The solar neutrino data includes a total of 118 



data points: the Galliu m |30l . l3lj and Chlorine [29(] (1 data point each) radiochemical rates 



the Super-Kamiokande [32| zenith spectrum (44 bins), and SNO data reported for phase 1 
and phase 2. The SNO data used consists of the total day-night spectrum measured in the 
pure D 2 (SNO-I) phase (34 data points) [33|, P ms the full data set corresponding to the 
Salt Phase (SNO-II) (34|. This last one includes the neutral current and elastic scattering 
event rates during the day and during the night (4 data points), and the charged current 
day- night spectral data (34 data points). It is done by a x 2 analysis using the experimen- 
tal systematic and statistical uncertainties and their correlations presented in [HI, together 
with the theoretical uncertainties. In combining with the SNO-I data, only the theoreti- 
cal uncertainties are assumed to be correlated between the two phases. The experimental 
systematics errors are considered to be uncorrelated between both phases. 

In the analysis of KamLAND, we neglect the effect of the long-range forces due to the 
small Earth-crust density in the evolution of the reactor antineutrinos. For KamLAND we 
have used the observed events as a function of L /E u , with L fixed at 180 km 3|| and 
minimized the x 2 function 
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10 3 10 2 10 1 1 10 

Au (TeV) 

FIG. 3: Bounds from the analysis of solar and KamLAND data on the fundamental parameters 
(A u , M) for a scalar unparticle for d = 1.1, 1.3, 1.6 and duv = 3. The region below the curves is 
excluded. We have fixed as = 10~ 3 . The dashed lines correspond to the approximation given in 

Eq. A25D- 



Xkl 



24 

£ 

i=l 



2{fRf -RT + K? log 



fR\ 



th 



l-f 

0.041 



(24) 



with respect to Am^, tan and /. Here Rf- and R\ are, respectively, the observed and 
theoretical (which depends of the oscillation parameters) number of events, for the i-th bin. 

One must also notice that the lagrangian given in Eq. ([3]) also leads to extra neutral 
current contributions to neutrino interaction cross-sections at the detector. However, we 
will ignore those effects which in the absence of cancellations will add-up to the propagation 
effect. In this sense the bounds established here are conservative. 

As an illustration we show in Fig. [3] the bounds obtained from our combined analysis of 
solar and KamLAND data on the fundamental parameters for an scalar unparticle. For the 
sake of concretness we have chosen as > but similar order of magnitude bounds can be 
derived for negative as. 

First we notice that the bounds from the long-range force effects are orders of magnitude 
stronger than the corresponding ones from the unparticle-induced invisible solar neutrino 



11] . This justifies our approximation of neglecting the decay in the 



decay described in Ref. 
present work. 

The order of magnitude of the bounds derived in Fig. [3] can be quantitatively understood 
from the comparison shown in Fig. [TJ between the unparticle induced potential function and 
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that expected from a standard finite-range interaction. As seen in the figure, for any value 
of d < 2 the W function is always larger than the corresponding one for a d = 1 force of 
range A = 0.1 R Q . Consequently, for any value of d < 2 the bounds on the unparticle effects 
must be comparable or stronger than the limits obtained in Ref. j26| for a leptonic scalar, 
vector, or tensor, long-range force with coupling ks, k>v an d kt, respectively, and of range 
A = 0.1R & with the identification 

, A \2(<fcv-i) 2 r [ d+ i] r [ d _i] 1 

K s -> «S TT -t 7T J ^4 ^TT^T < 6.8 x 10~ 45 , 25 

s \m) 7T r[2d] (i^Aj 2 ^- 1 ) - 1 ; 

« fM 2(<w_1) 2 E^+llEEzi! 1 < 45 xio- 53 

mp (K\ 2{duv - 1] 2 r[d + i]r[rf-l] 1 2 x 10-61 

" T ~* ^-i T[2d] (R Q A U ) 2 ^ ~ eV ' {Z<} 



For comparison we also plot in Fig. [3] the bounds directly given by Eq. (j25p . As seen the re- 
sults from the full analysis are in reasonable quantitative agreement with the approximations 
above. 

From Eqs. (1251) - (J27j) we can see that for a given value of A u at the same scale dimension 
d, the bounds on M are ~ 100 times more stringent for vector unparticles. For tensor 
unparticles the bounds on M will be ~ 3/ ^/A„/TeV times those for the scalar case. The 
bounds tend to get substantially relaxed as d increases. Consequently, , if the low energy 
scale is of the other A u ~ TeV, one cannot provide very significant constraints on vector and 
tensor unparticles arising in a full conformal symmetry - for which d > 3, 4 respectively- 
from their long-range effects on neutrino oscillation data. 



V. DISCUSSION 

We have studied the effect of the long-range forces induced by the coupling of unparticle 
operators to fermions of the SM, which are exerted along the the neutrino trajectory from 
its production point inside the Sun to its detection at the Earth. We have used data from 
solar neutrino experiments and the newest KamLAND results to show that one can place 
stringent constraints on the parameters involved in the description of these new interactions. 

For a scalar unparticle, for instance, assuming the non-universality of the lepton couplings 
to unparticles to be of the order of a few per mil, i.e. a s = 10~ 3 , and c?uv = 3, we find 
that the mass M of the UV exchanged particle is constrained to be M > 8 x 10 8 TeV 
(M > 7 x 10 9 TeV), for the scaling dimension d = 1.1, if the lower energy scale is at A u = 
1 TeV (10 TeV). 

In general the bounds derived here are much stronger than the reach at present and near 
future colliders. They are comparable with those imposed from cosmology and astrophysics 
and from the effect on the modification of the Newtonian law of gravity. For the sake of 
comparison we show in Table [I] our bounds together with those imposed by BBN [l3[ , and 



from the compilation in Table 2 of Ref. [15j of the bounds from astrophysics and 5th force 
experiments. In order to make this comparison we have assumed all Cg to be of order 1 but 
that some violation of universality is allowed at the per mil level. 

We can see that for low values of the scale d, the bounds derived in this work are com- 
parable to the most stringent ones, coming from scalar unparticle mediated 5th-force tests. 
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Scalar Unparticle 


M > (TeV) 




1 TeV 


10 TeV 


d 


1.3 


1.6 


1.3 


1.6 


This work 


1.3 x 10 6 


100 


9 x 10 6 


500 


BBN 


3.2 x 10 5 


5.8 x 10 3 


4.7 x 10 6 


1.2 x 10 6 


Eotvos-type 


1.9 x 10 6 


3.5 x 10 3 


1.3 x 10 7 


1.6 x 10 4 


Energy loss from stars 


800 


26 


5.5 x 10 3 


120 


SN 1987A 


430 


55 


3 x 10 3 


260 



TABLE I: Limits on M from various sources testing for signals from scalar unparticle couplings to 
fermions. Here d\jy = 3 was used in deriving all bounds and a s = 10~ 3 was used for our work. 

However, unlike these last ones, the bounds from the solar neutrino oscillation effects would 
apply even if the unparticles do not couple baryons. 

Conversely, the bounds we derived in this work can also be converted into severe con- 
straints on universality violation of the neutrino couplings to unparticle operators with 
scales which can be accessible at future colliders. For example, for M = A u = 1 TeV and 
C f s = 0(1), the bounds in Eq. Q25D imply that 

C ^~^ a < 4.5 x lCr 38 - 1.5 x 10" 9 (28) 

for d = 1.1—1.6. Therefore in the absence of a model justification of why flavour effects 
on the neutrino couplings to unparticle operators are more constrained than on unparticle 
couplings for other fermions, this analysis implies that most likely the fermion couplings to 
unparticles which can be tested at colliders will be flavour blind. 
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